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Abstract. Many low- mass stars in the Orion nebula are associated with 
very compact (~ 1 arcsec) emission knots, known variously as proplyds, 
PIGs or LV knots. Some of these knots are teardrop-shaped, with "tails" 
pointing away from the massive star 6^ Ori C, which is the principal ex- 
citing star of the nebula. We discuss models of such knots, which invoke 
the interaction of the fast stellar wind from 9^ Ori C with a transonic pho- 
toevaporated flow from the surface of an accretion disk around a young 
low-mass star. We review previous analytic work and compare the results 
of the model with the observed brightnesses, morphologies and emission 
line profiles of the knots, as well as presenting new results from numerical 
hydrodynamical simulations. 



1. Introduction 

The proplyds are bright compact emission line knots, with sizes of order 
0.5-2.0 arcseconds, that are found in the inner region of the Orion nebula 
(Laques and Vidal, 1979; Garay et al., 1987; Churchwell et al., 1987; Felh 
et al., 1993; O'Dell et al., 1993; O'Dell and Wen, 1994; McCaughrean, 1997) 
and nearly all of which contain an embedded low-mass star (Meaburn, 1988; 
McCaughrean and Stauffer, 1994). Many of the proplyds show a head/tail 
morphology, in which the tail points away from the star 6^ Ori C, the most 
massive star of the Trapezium cluster. Emission line spectroscopy of the 
proplyds in the [O iii] 5007A line (Massey and Meaburn, 1993; Massey and 
Meaburn, 1995; Henney et al., 1997) show a bright central core with full 
width half maximum (FWHM) of ~ 50 km s~^, together with faint wings 
extending out to > 100 km s~^ from the line center. 
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The obvious explanation for these objects is that the radiation and stel- 
lar wind from 9^ Ori C is interacting with the circumstellar material around 
the low-mass stars. In our models we assume that this circumstellar mate- 
rial is in the form of an optically thick, geometrically thin accretion disk. 
The effect of the radiation from 9^ Ori C will be to ionize the material in 
the disk, which will produce a photoevaporated flow away from the disk 
and towards the ionizing source. Hence, to produce the tails pointing away 
from 9^ Ori C, this flow must somehow be confined and redirected. Two 
possible candidates for this confinement mechanism are the exciting star's 
radiation pressure and the ram pressure of its stellar wind. However, only 
the second of these is feasible, as is shown in § 2. In § 3 a simple analytic 
model of the resultant two-wind interaction is briefiy described and the 
successes and failures of the model in explaining observed properties of the 
proplyds are outlined in § 4. In § 5, preliminary results of numerical hy- 
drodynamic simulations are presented, which remove some of the arbitrary 
assumptions of the model. Complications such as the possible existence of 
a neutral photodissociated flow are critically discussed in § 6. 

2. Confinement Mechanisms — Radiation vs. Ram Pressure 

The gas pressure at the base of the ionized flow can be calculated simply 
by equating the numerical flux, Fq, of Lyman continuum (Ly-c) photons 
arriving at the ionization front (IF) with the numerical flux of newly-ionized 
ions entering the photoevaporated wind: 



where no is the ion number density and uq the ion velocity at the base of 
the wind. This leads to the following expression for the gas pressure 



where fi is the mean atomic mass (~ 1.3), mn is the mass of hydrogen, cq 
is the sound speed in the ionized gas (~ 12 km s""*^) and is the Mach 
number at the base of the flow, which will be of order 1-2 (Dyson, 1968; 
Kahn, 1969; Bertoldi, 1989). 

The unattenuated ionizing flux is given by F^, = S^/4Tr(f, where S"* is the 
stellar ionizing photon rate, for which estimates vary between 7 x 10^^ 
(Panagia, 1973) and 3 x lO'^^s"^ (Bertoldi and Draine, 1996), and d is the 
distance of the proplyd from the exciting star. However, at the distances of 
the proplyds from 9^ Ori C, most of this flux is used up in maintaining the 
ionization state of the photoevaporated flow against recombination. With 
the assumption that Fo <^ F^, one can write 



Fq = noUo , 



(1) 



(2) 
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where tq is the dust absorption optical depth of the flow, is the Case 
B recombination coefficient (2.6 x 10~^^ cm^ s~^), ra is the disk radius 
and /(tq) — (3 + tq)~^ depends slightly on the assumed geometry (Hcnney 
et al., 1996). Assuming tq = 0, one finds (Eqs. 1 and 3) that the percentage 
of ionizing photons reaching the ionization front is 

/3% = ^ = 1.5Modlrr,-^/^^4V/^ (4) 

where dn and ri5 are the proplyd distance and disk radius measured in 
units of 10-*^^ cm (~ 0.3 pc) and 10-*^^ cm (~ 66 au) respectively and 549 is 
the stellar ionizing photon rate in units of 10^^ s~^. Allowing for the effects 
of dust makes little difference to this estimate. 

The ionizing radiation pressure from 6^ Ori C that acts on the photoe- 
vaporated flow can be written as 

P,. = ^. (5) 

where c is the speed of light, is the unattenuated ionizing flux from 
9^ Ori C, and {hu) is the mean energy of ionizing photons absorbed in the 
flow (~ 13.6eV). Hence, 

" ^i;;^^^^ = 0.033 A^o/3^„ = 0.022 d,, r^, S^, . (6) 

This ratio is always significantly less than unity, therefore the ionizing radi- 
ation pressure is incapable of confining the photoevaporated flow. If there 
were enough dust opacity at the base of the flow, then it is conceivable 
that the non-ionizing radiation from 6^ Ori C (FUV, optical) may make a 
significant contribution to the radiation pressure. However, the bolometric 
luminosity of Ori C is only ^ 3 times its Ly-c luminosity, so the above 
conclusion is unchanged and radiation pressure falls an order of magnitude 
short of the thermal pressure even for the closest proplyds {du ~ 0.5). 

Turning now to the ram pressure, -Phyd) of the stellar wind from 9^ Ori C, 
this will be given by 

Phyd = Pwvl , (7) 

where and arc respectively the stellar wind density and velocity. 
Since the wind is radiation-driven, one would expect the ratio -Phyd/^'rad to 
be of order unity and, using the observed parameters of 9^ Ori C (Howarth 
and Prinja, 1989; Panagia, 1973), one finds (Henney et al., 1996) that this 
is indeed the case. However, although the radiation pressure must act on 
the base of the wind where the majority of the recombinations occur, the 
ram pressure need not do so, but will act on the surface of contact between 
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Figure 1. Schematic diagram of the two- wind interaction model. Ionizing photons from 
6^ Ori C drive a photoevaporated wind from one face of the circumstellar disk, which 
interacts supersonically with Ori C's stellar wind. 



the evaporated flow and the stellar wind, wherever that may be. Hence, 
the photoevaporated gas will flow divergently away from the disk until its 
pressure (reduced by geometric dilution) falls to that of Phyd; at which 
point it can be confined by the stellar wind. 

3. Analytic Two- Wind Models 

The analytic model (Henney et al., 1996) depends chiefly on the dimension- 
less parameter A = -Pgas/-fhyd- From the discussion of the previous section, 
it is evident that A > 1, in which case the photoevaporated flow, which is 
initially mildly supersonic (Dyson, 1968; Kahn, 1969; Bertoldi, 1989), will 
begin to flow freely away from the disk. It is assumed that the streamlines 
are straight and that the initial flow diverges with a half-opening angle of 
45°. If the velocity remained constant, the density would fall as (1 -|- -z)"^, 
where z is the height above the disk in units of the disk radius, but a 
pressure gradient causes the flow to accelerate. 

The flow will shock at the point where its pressure has fallen to that of 
-Phydj which occurs at a distance 



D ~ 



1.19 (In A) 



V4;^i/2 



rd = 4-20 rd , 



(8) 
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where Oq is the inchnation angle of the disk normal with respect to the 
direction to 9^ Ori C and, for the second equality, A = 10-100 is assumed. 
The shock will be radiative, so can be treated as isothermal, although the 
radiation from the shock itself makes a negligible contribution to the pro- 
plyd luminosity. A shock will also form in the wind from 9^ Ori C, but this 
will be non-radiative, hence the assumption of ram pressure balance used 
to derive equation 8 is not strictly valid (see § 5). 

A thin, almost flat, layer (Mach disk) of shocked photoevaporated flow 
material forms parallel to the circumstellar disk and gas flows outwards 
along this layer, reaching a velocity at the edge of 

ve ~ 3.8(ln A)^/2 Co = 70-100 km s'^ . (9) 

The gas is then swept back by the wind of Ori C to form the pro- 
plyd wings and tail. Figure 1 illustrates the components of the model in a 
schematic form. 

For reasonable values of A, the photoevaporated disk wind is the bright- 
est component of the model (with a luminosity ~ 0.5A^/^ times that of the 
Mach disk plus tail) and also the smallest, leading to a core-halo morphol- 
ogy (Henney et al., 1996, Fig. 11). 

4. Comparison with Observations 

The ensemble properties of the proplyds are quite well reproduced by the 
analytic model. The models show good agreement with the observed trends 
of proplyd size and luminosity vs. distance from 9^ Ori C, the former in- 
creasing and the latter decreasing (McCullough et al., 1995; Henney ct al., 
1996, Figs. 9 and 10). The implied circumstellar disk radii are between 20 
and 60 au (ris ~ 0.3-1.0). These correlations, however, are rather insensi- 
tive to the details of the model. 

The morphologies of individual proplyds are compared with model pre- 
dictions in Figure 2 (a larger sample is given in Fig. 14 of Henney et al., 
1996), where it can be seen that the models successfully reproduce sin- 
gle and double tails, both of which are observed (O'Dell and Wen, 1994; 
O'Dell and Wong, 1996; Johnstone et al., 1996). However, the double tails 
may be merely the result of absorption in the core of the tail, which is 
not consistent with the models as they stand. The crescent head observed 
in many proplyds would correspond to the Mach disk in the models, but 
this is rather problematic since the models predict that this should be less 
bright than the photoevaporated wind component (§ 3) , which is not the 
case for most proplyds, although dust absorption at the base of the wind 
(r ~ 0.5-1 for the closest proplyds) would alleviate this problem (this is 
included in the fit to OW 158-327). 
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Detailed comparisons between model predictions and high resolution 
[O III] 5007A spectra of individual proplyds are presented in Henney et al. 
(1997). The evaporated wind produces the bright core of the line, with 
width of a few times the sound speed in the ionized gas, while the Mach disk 
and tail produce the high-velocity line wings that are observed. However, in 
order to reproduce the ~ 100 km s^^ widths of the line wings seen in LV 5 
(OW 158-323) and LV 2 (OW 167-317), values of A = 50-200 are required, 
which are 3-4 times larger than those found in fitting the morphologies of 
the same objects (Henney et al., 1996). 

5. Hydrodynamical Simulations 

Figure 3 shows the results of an example numerical simulation of the two- 
wind interaction (Henney and Arthur, 1997). In this simulation, the cir- 
cumstellar disk (oriented vertically in the figure) is inclined by 45° with 
respect to the direction of the stellar wind from 9^ Ori C (other parameters 
are described in the figure caption). The transfer of ionizing radiation in 
the photoevaporated flow is not calculated sclf-consistcntly in the models, 
but the boundary conditions are assumed constant over the disk surface 
and are taken from the analytic model. Also, the simulation parameters 
correspond to a rather small value of A (~ 8), since a larger value would 
require an unfeasibly large computational grid. 

The main differences with respect to the analytic calculation are due to 
the relaxation of two arbitrary assumptions of the model. Firstly, the pho- 
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Figure 3. Hydrodynamical simulation of the two-wind interaction, calculated in 2D 
slab symmetry with a grid size of 300 x 300 cells. Grayscale shows the gas density while 
contours show the gas pressure (both logarithmic scale). Arrows show gas velocity. Pho- 
toevaporated disk material (white arrows) has an isothermal equation of state. Stellar 
wind material (black arrows) has an adiabatic equation of state. Note the different ve- 
locity normalizations of the black and white arrow lengths. The parameters of the model 



are ra = 30 au, 6o = 4:5° , no = 2 x 10 cm nw = 1cm 



= 1000 km s-^ 



toevaporated flow is calculated self-consistently, instead of being assumed 
to follow straight streamlines with opening angle 45°. Pressure gradients 
in the mildly supersonic wind in fact cause the flow to diverge increasingly 
with distance from the disk. Secondly, the (non-radiative) shock in the 
stellar wind is treated properly, instead of merely assuming ram pressure 
balance between the two winds as was done in the analytic model. Both 
these factors affect the shape of the dense layer of shocked photoevaporated 
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wind material, which is more curved than in the analytic calculation. The 
flat "Mach disk" of the analytic model is no longer apparent and there is no 
sharp distinction between the Mach disk and tail. The gas velocities reached 
in the shocked photoevaporated wind are also slightly smaller than in the 
analytic model. Unfortunately, the use of slab symmetry, which allows the 
asymmetric interaction to be modelled in two dimensions, means that it 
is not possible to produce emission maps or spectra from the simulations. 
Nonetheless, the morphology and kinematics of the simulations are rather 
similar to those of the analytic calculation modulo the differences noted 
above. The arguments for and against the two- wind interaction model are 
hence little affected. 

6. Discussion and Speculation 

Despite the success of the two-wind interaction models in reproducing the 
observed morphologies and kinematics of the proplyds, some problems re- 
main. In particular, the apparent absorption in the tails of some objects 
and the [Oni]/IR arcs (Bally et al, 1995; Hayward et al., 1994) that 
are seen between the closer proplyds and 9^ Ori C are both hard to ex- 
plain with the two- wind model. An alternative view (Bally et al., 1995; 
Johnstone et al., 1996) is that the disk evaporation is controlled by non- 
ionizing FUV photons, with the ionization front occurring away from the 
disk. The proplyd morphology would then be determined by the shape of 
the ionization front. The hydrodynamic interaction with the stellar wind 
from 6^ Ori C would, on this view, still occur, but farther out in the flow, 
perhaps producing the [Oiii]/IR arcs. This model has had most success 
in explaining the characteristics of HST 10 (OW 182-413), but this object 
does not seem to be typical of proplyds as a class (in particular, its tail does 
not point exactly away from 9^ Ori C and it may not contain a central star). 
Johnstone et al. (1996) compare the mass-loss rates from photodissociated 
and photoionized disk winds and conclude that the former will dominate for 
all proplyds. However, they assume that the warm (~ 1000 K) photodisso- 
ciated gas will be able to flow freely away from the disk at its sound speed 
(~ 2.5km s~-^), but this is not necessarily the case. 

If one allows, for the sake of argument, that a free-flowing photodissoci- 
ated wind, with a density at its base of 10®nn,6 cm~^, is initially established, 
then, once the ionizing radiation from 9^ Ori C is switched on, an R-type 
ionization front will be driven rapidly into the flow. For proplyds closer 
than 

d;7^0.55n->r5'^'54f , (10) 

the flow will be immediately ionized all the way down to the disk. For 
proplyds further away, the ionization front undergoes a transition to D- 
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type some way out from the disk, at which point its progress wih slow and 
it will begin to drive a shock into the atomic flow. The density of newly 
ionized gas rig will adjust itself to that given by equations 1 and 4 of § 2, 
but with Td replaced by the radius of the ionization front. The pressure of 
the ionized gas will be roughly 20 times that of neutral gas of the same 
density, so that for proplyds closer than 

d'(,c,20n-},r^,'/'sli\ (11) 

the shock will reach the surface of the disk before stalling, hence quenching 
the neutral flow in a time of 30-500 years. In proplyds farther away than 
(i'/y from 6^ Ori C, the shock will stall at a distance zor^ from the disk, 
where 

1 + 0.136 nl(l dli' rif S^.'^^ . (12) 

The chief uncertainty in the estimates of d'^ and d'/^ is the density at 
the base of the neutral flow. However, taking the parameters of HST 10 
(Johnstone et al., 1996), ris = 1.3, zq = 2.3, dn = 5, and assuming nn,6 is 
the same for all proplyds, one finds that d[j ~ 0.03 rj"5^^^ and dij ~ r^^^"^- 
No proplyds are observed with dn < d'^ but a substantial fraction have 
dn < d'{-j, although the exact number depends on the distribution of disk 
radii. This can only be determined directly for the dark silhouette disks 
(McCaughrean and O'Dell, 1996), which show ris = 0.4-7.6, but the bright 
proplyds are likely to have smaller disks (ris ~ 0.1-1.0, Henney et al., 1996; 
Johnstone et al., 1996). Hence, roughly half of all bright proplyds will not 
have an extended neutral evaporated flow. 

The real situation is undoubtedly much more complicated than por- 
trayed above (c.f. Bertoldi and Draine, 1996), but the basic argument, that 
the neutral flow must have a higher pressure than the overlying ionized flow 
in order to exist, should remain valid. A further problem for the neutral 
flow is gravity: the escape speed from the circumstcUar disk will equal the 
sound speed in the neutral gas at a disk radius of resc,i5 — 2.1M*, where 
M^, is the mass of the central star in solar masses (~ 0.1-2, McCaughrean 
and Stauffcr, 1994). Hence, except for the proplyds with the lowest mass 
central stars, gravity will dominate the dynamics of the photodissociated 
region. 

Note that the argument against radiation pressure in § 2 applies a for- 
tiori to a neutral photodissociated flow since its pressure would have to be 
larger than the ionized flow. However, the confinement problem could be 
circumvented if it were maintained that the material in the tail, instead of 
having been redirected from an initial flow towards 9^ Ori C, was instead 
part of a flow from the back side of the disk, possibly driven by the diffuse 
radiation field. This could also explain the absorption seen in the core of 
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some tails, but whether the flow would be dense enough for this is not clear. 
Alternatively, the tails could be formed from the remnants of a dense slow 
wind from a massive star (Sutherland ct al., 1997). 

In conclusion, the two-wind interaction model has had qualified success 
in explaining the observed properties of the proplyds closer to 6^ Ori C. 
Various discrepancies remain, however, and further work is necessary both 
in extending this model and in developing alternatives. 
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